Ground state of a double-exchange system containing impurities: bounds of 
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We study the boundary between ferromagnetic and non-ferromagnetic ground state of a double- 
exchange system with quenched disorder for arbitrary relation between Hund exchange coupling and 
electron band width. The boundary is found both from the solution of the Dynamical Mean Field 
Approximation equations and from the comparison of the energies of the saturated ferromagnetic 
and paramagnetic states. Both methods give very similar results. To explain the disappearance 
of ferromagnetism in part of the parameter space we derive from the double-exchange Hamiltonian 
with classical localized spins in the limit of large but finite Hund exchange coupling the t — J model 
(with classical localized spins). 
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The problem of ferromagnetism in the double-exchange 
(DE) model |]-|, which appears due to the Hund ex- 
change coupling between the core spins and the mobile 
carriers, has a long history (see Ref. Q| and references 
therein). Frequently the modeling deals with the case 
of infinite Hund exchange. In this approximation the 
ground state of the system is always ferromagnetic (FM) . 
To describe possible non-FM ground states of such a sys- 
tem one has to introduce direct antiferromagnetic ex- 
change between the core spins. However, the finite Hund 
exchange by itself generates effective antiferromagnetic 
exchange, thus destroying ferromagnetism in a part of 
the system parameter space. 

In our previous publication [^j on the basis of DMFA 
we obtained the closed formula for the ferromagnet- 
paramagnet (FM-PM) transition temperature T c for the 
double-exchange system for arbitrary relation between 
Hund exchange coupling and electron band width. In 
this paper we present a detailed study of of the bound- 
ary between FM and non-FM ground state of the system 
found from that formula, and, independently, from com- 
parison of energies of the saturated FM and PM states. 
In addition, from the double-exchange Hamiltonian with 
classical localized spins in the limit of large but finite 
Hund exchange coupling we obtain the classical version 
of the t — J model. 



II. HAMILTONIAN AND DMFA EQUATIONS 



Consider the DE model with random on-site energies. 
The Hamiltonian of the model is 



_ J ^2 m ™ ' a af3 C ia c nl3, (1) 
na/3 

where t n —n' is the electron hopping, J is the effective 
exchange coupling between a core spin and a conduc- 
tion electron, a is the vector of the Pauli matrices, and 
a, (3 are spin indices. We express the localized (classical) 
spin by m„ = (m n x ,m n v ,m n z ) with the normalization 
|m| 2 = 1. To take into account the chemical disorder 
introduced by doping impurities, which is generic for the 
manganites and many other DE systems, we consider the 
case random on-site energy V n . 

In a single electron representation the Hamiltonian can 
be presented as 



(V n - Jm n ■ a) Sr, 



(2) 



the first is translationaly invariant, the second describes 
quenched disorder, and the third - annealed disorder. 

The DMFA, as applied to the problem under consider- 
ation, is based on two assumptions. The first assumption 
is that the averaged, with respect to random orientation 
of localized spins and random on-site energy V, locator 



Gioc(z) = (G nn {z) 



m,V 



where 



G(z) = (z-HY 



(3) 



(4) 



can be expressed through the the local self-energy S by 
the equation 
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ioc(z) =9o(z- . 
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where 
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is the bare (in the absence of the disorder and exchange 
interaction) locator. Thus introduced self-energy satisfies 
equation 



(jm ■ a - V n + t(z s ) 



(13) 



G\ oc (z) 



1 



G^ c (z)+t(z)-V n + Jm-a i 



(7) 
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The system of equations (|[) and (0) is very much similar 
to the well known CPA equations (see || and references 
therein), as generalized to the case when the quantities 
G,T, and g are 2x2 matrices in spin space 0. The 
system of equations however, is not yet closed. The av- 
eraging with respect to annealed disorder is principally 
different from the averaging with respect to quenched 
disorder. 

The second assumption of the DMFA is the prescrip- 
tion for the determining, in our case, the probability of 
a spin configuration self-consistently with the solutions 
of the Eqs. (|) and (@). To formulate the DMFA equa- 
tion for this probability, taking into account both kinds 
of the disorder, let us start from the general formula for 
the partition function 

Z Vn = J exp ^-Tr^logG^ \[dm n , (8) 

where z s = ioj s + /i; lo s is the Matsubara frequency and 
H is the chemical potential. The averaging over {m n } is 
given by 

($)m " ^;/ exp (- Tr E lo s^)) $ HlI rfm - ( 9 ) 

All observables, in particular thermodynamic potential 
f2, should additionally be averaged over the realizations 
of the quenched disorder; in particular 



is the change of the thermodynamic potential of the elec- 
tron gas described by the Green's function G\ oc due to 
interaction with a single impurity [p|jic|l . 

The right hand side of Eq. (|l2|), is a complicated non- 
linear functional of Py(vsx). However, if we are interested 
only in the transition temperature T c , the problem can 
be reduced to a traditional mean field (MF) equation. In 
linear with respect to magnetization M approximation 
Eq. (|l^) takes the form 



P Vn (m) oc exp (-/3I Vn M • m) . 
Non-trivial solution of the MF equation 



(14) 



M = J (P Vn (m)) v mdm. (15) 
can exist only for T < T c , where T c = i (Iv„)v 

III. T c FOR THE SEMI-CIRCULAR DOS 

For simplicity consider the semi-circular (SC) bare den- 
sity of states (DOS) N (e), the bandwidth being 2W. 
Then 



9o(z) = 



N (e)de _ 2 
z-e ~ W 
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For this case 



S = z — 2wG\ oc — G 



loc' 



(16) 



(17) 



where w — W 2 /8. Thus from Eqs. and (Q) we obtain 
a single equation for Gi oc 



(10) 



Gloc(z) 



z — 2wG\oc(z) — V n + Jm ■ o 



(18) 



The DMFA approximates the multi-spin probability 
Zy 1 exp (^— TrlogG^j as a product of one-site probabil- 
ities in such a way, that 

sn , , 

,0. (11) 



SG 



loc 



The result for the one-site probability reads (for details 
of the calculation see Ref. S): 



P Vn (m) oc exp [-/?AQ m) yJ 



(12) 



where 



and Eq. (|T^) can be presented as 

Att m y n = 4 E lo gdet z s - 2wG loc (z s ) 

s 

-V n + Jm-a]e lUsG +. (19) 

In linear with respect to M approximation 

Gioc = gl - hJM ■ a, (20) 

where g is locator in paramagnetic phase, given by the 
equation 
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z — 2wg 
1 



V n -J 



z — 2wg — V n + J 
and the quantity h is given by the formula 

<A 2 V \ -2 w <A y J y ' 



4J 2 M) 



where 



Ay„(z s ) 



1 



[z. - 2wg(z s ) - V n ] - J 2 



(21) 



(22) 



(23) 



Expanding Eq. (|19| ) we obtain the effective exchange in- 
tegral is 



AJ 2 w 



^h(z s )A Vn (z s ). 



(24) 



If we transform the sum over the imaginary Matsubara 
frequencies in the right-hand side of Eq. ( p4| ) to integral 
over real energies E, we obtain for the T c 



AJ z w 
3tt 



f(E)Im[h(E + )(A Vn (E + )) v ]dE, (25) 



where f(E) is the Fermi function, and E + = E + iO. 
To the best of our knowledge this closed formula which 
takes into account both finite value of Hund exchange and 
quenched disorder was obtained and studied by us for the 
first time j^] . According to this formula the boundary be- 
tween the FM and non-FM ground states is found from 
equating T c to zero, which reads 



lm[h(E+) (A Vn (E+)) v }dE = 0; 



(v) 

the Fermi energy Ep is found from the equation 



(26) 



Im g(E+) dE, 



(27) 



where n is the number of electrons per site. 

Thus found T c is the temperature at which the para- 
magnetic state became unstable with respect to small 
spontaneous magnetic moment. This becomes especially 
obvious, if we obtain the equation for the T c by con- 
structing Landau functional || . Independently, one may 
consider another approach to the same problem based on 
the comparison of the energies of the saturated FM state 
and PM state. The energy of the PM state is found from 
the equation 



E p = 

7T 



E Im g(E+) dE; 



(28) 



the energy of the saturated FM state and the appropriate 
Fermi energy are found from the equations: 

E (f) 

E f = -- I " E [Im 91 (E + ) + Im 9l (E + )] dE 

E (f) 

n = -- [ ' Im [Im 9] (E + ) + Im gi (E+)}dE, (29) 

7T J-oo 



where the FM locator is given by the equation. 

! \ . 

, z - 2wg -VnT J I v 
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(30) 



The boundary of the ferromagnetic region at the phase 
diagram can be found from the equality between the two 
energies: E p — Ef. We thus use an alternative to the 
second assumption of the DMFA. 



IV. FM - NON-FM BOUNDARY 

The boundary of the FM ground state of the system 
in case of no quenched disorder is presented on Fig. 1. 
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FIG. 1. The FM - non-FM boundary for the case of no 
quenched disorder in the coordinates of relative strength of 
the Hund exchange J/W and electron concentration n. 

The diagram agrees with those obtained on the basis 
of numerical calculations [ pd] | and from qualitative rea- 
soning pQ |. 

To consider the influence of chemical disorder we con- 
sider the model in which V n = V with the probability 
x, and V n = with the probability 1 — x, thus x being 
the concentration of impurities. In this model we con- 
sider the number of electrons n and the concentration 
of impurities x as two independent parameters. Solving 
equation for the locator we obtain the boundaries which 
are presented on Figs. 2-4. It is interesting that ferro- 
magnetism is now precluded in much larger region of the 
J/W — n plane. 
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FIG. 2. The FM - non-FM boundary for V/W = .5 and 
x = .3 in the coordinates of relative strength of the Hund 
exchange J/W and electron concentration n. 
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FIG. 3. The FM - non-FM boundary for V/W = 1 and 
x — .3 in the coordinates of relative strength of the Hund 
exchange J/W and electron concentration n. 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

n 

FIG. 4. The FM - non-FM boundary for V/W = 2 and 
x — .3 in the coordinates of relative strength of the Hund 
exchange J/W and electron concentration n. 



The boundaries between the phases make sharp bends 
for static disorder strong enough to create a gap between 
conduction and impurity bands both in ferromagnetic 
and paramagnetic phases . In this case the number of 
states in the impurity band is equal to the concentration 
of impurities, and for n = 1 — x the Fermi energy is in the 
gap. To emphasize this fact we plot the density of states 
both for FM and PM phases at the values of parameters 
V/W, x and J/W which correspond to the sharp bends 
on Figs. 3 and 4. 
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FIG. 5. Density of states for V/W = l,x= .3, J/W = 1.53 
for saturated FM phase (left) and PM phase (right); occupied 
states for n = 0.7 are filed. 
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FIG. 6. Density of states for V/W = 2,x= .3, J/W = 2.53 
for saturated FM phase (left) and PM phase (right); occupied 
states for n = 0.7 are filed. 



APPENDIX A: FROM THE 
DOUBLE-EXCHANGE HAMILTONIAN TO THE 
T — J MODEL: CLASSICAL SPINS 

The Hamiltonian of the DE model is 



1 



H 



Jh m„ ■ (J a ^ na , c n& 
na(3 



nn' a 



(Al) 



where t n - n i is the electron hopping, Jh is the effective 
Hund exchange coupling between a core spin and a con- 
duction electron, a is the vector of the Pauli matrices, 



4 



and a, j3 are spin indices. We express the localized (classi- 
cal) spin by a unit vector whose orientation is determined 
by polar angle and azimuthal angle cf). In a single elec- 
tron representation the Hamiltonian can be presented as 



H nn ' = H e 



H 



— J/fm n • a5 nn ' — t, 



(A2) 



We consider the case of strong exchange: Jh ^ W where 
W is the electron band width. In this case we should first 
diagonalize the exchange part of the Hamiltonian. This 
is done by choosing local spin quantization axis on each 
site in the direction of m. In this representation the 
Hamiltonian is 11131 



Hi, 



where 



°n'n a r< 



(A3) 



"n "n 1 . . "n . "ri UA, ,—(h„) 

a-nn' — cos — cos h sin — sin — e w Vn> 

nn 2 2 2 2 



Kn> = sin — cos — e 



*t>n _co S ^i sin Z2l e -^. (A4) 
2 2 x ' 



The transformation to local spin quantization axis in- 
cluding an additional Euler rotation angle, which leads 
to a more involved effective Hamiltonian than A3, was 
introduced by Nagaev . 

The next step, like it is done in the derivation of the 
t — J model from the Hubbard model , is to apply a 
canonical transformation 

H ->H = e s He- s = H+[S,H} + -[S[S, H]] + ... (A5) 

which excludes all band-to-band transitions. This can be 
achieved if we chose the operator S in the form 



tn 



We have 



[S, [S,H° 



2J H 



[S, H e 



7 * 

"Vn a n'n 



t, 



(A6) 



-\bnn" 





b nn > 

K'n 

1 

\b nn "\ 



(A7) 



where J nn " = \t n n" \ 2 /{2Jh)- Keeping terms up to the 
second order with respect to t (and only site-diagonal 
part of the second order terms) we obtain 



nn' L nn 



Inn' 
Qn'n 

^2 Jnn"(m n ■ m„" - l)a z 5 

n't 



(A8) 



In the second quantization form the Hamiltonian (AS) 
has the form (ignoring the constant term) 



nn' 

Jnn'TXln ' m„/ (1 - d n d n ), 



(A9) 



where for the case of hole doping dJ(d) is the operator 
of creation (annihilation) of the hole, and for the case of 
electron doping it is the operator of creation (annihila- 
tion) of the electron. 

Looking at the Hamiltonian obtained, we see that large 
but finite Hund exchange dynamically generates anti- 
ferromagnetic exchange (the second term in Eq. flAS|). 
This term, however, is not independent upon the elec- 
tron (hole) subsystem. In this regard this Heisenberg like 
term resembles the first non-Heisenbergian term (kinetic 
exchange) . 
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